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Abstract. Vacuum bubbles nucleate at rest with a certain critical size and subsequently 
expand. But what selects the rest frame of nucleation? This question has been recently 
addressed in [1] in the context of Schwinger pair production in 1+1 dimensions, by using a 

4Ih model detector in order to probe the nucleated pairs. The analysis in [1] showed that, for 

a constant external electric field, the adiabatic "in" vacuum of charged particles is Lorentz 

^-H invariant, and in this case pairs tend to nucleate preferentially at rest with respect to the 

^ detector. Here, we sharpen this picture by showing that the typical relative velocity between 

^"^ —1/3 

00 the frame of nucleation and that of the detector is at most of order Av ^ S^ ' <C 1. Here, 

)^ 5^ ^ 1 is the action of the instanton describing pair creation. The bound Av coincides with 

the minimum uncertainty in the velocity of a non-relativistic charged particle embedded in 
a constant electric field. A velocity of order Av is reached after a time interval of order 

—1/3 

CO At ^ 5^ ro <C ro past the turning point in the semiclassical trajectory, where rg is the 

^ size of the instanton. If the interaction takes place in the vicinity of the turning point, the 

"^ semiclassical description of collision does not apply. Nonetheless, we find that even in this case 

S> there is still a strong asymmetry in the momentum transferred from the nucleated particles 

?-H to the detector, in the direction of expansion after the turning point. We conclude that the 

correlation between the rest frame of nucleation and that of the detector is exceedingly sharp. 
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1 Introduction 

False vacuum decay at zero temperature was first described in the context of field theory 
in a classic paper by Voloshin, Kobzarev and Okun [2]. These authors proposed that, as a 
result of quantum tunneling, true vacuum bubbles would nucleate at a certain rate per unit 
volume. The bubbles would nucleate at rest, with a critical size rg, and then expand with 
constant proper acceleration Tq" , due to the pressure difference between the false and the 
true vacuum. 

On the other hand, it was also pointed out in [2] that, due to the Lorentz invariance of 
the false vacuum, bubbles would not have any preferred reference frame in which to nucleate 
at rest. At first, this observation seemed to suggest that the total decay rate per unit 
volume should include an integral over the Lorentz group, in order to account for all possible 
rest frames of nucleation. Of course, such integral would be divergent, or at least cut-off 
dependent if a regulator was imposed. However, soon after this idea was proposed, Coleman 
[3] developed the instanton approach to false vacuum decay, showing that the decay rate is 
finite - and that integration over the Lorentz group does not play any role in calculating 
it ^. Coleman's argument, however, did not shed much light on a related and somewhat 
mysterious aspect of bubble nucleation, which requires clarification. If it is true that a long 
lived metastable false vacuum is approximately Lorentz invariant, what is it that determines 
the rest frame in which the critical bubble nucleates? 



^The reason is that instanton in this case is 0(4) invariant, and its analytic continuation (describing the 
bubble after nucleation) is invariant with respect to Lorentz boosts. Because of that, the final state in the 
asymptotic future is independent of the rest frame in which the critical bubble nucleates. Integrating over the 
Lorentz group would then amount to overcounting the final states. For a recent discussion of related issues, 

see [4-6]. 
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Figure 1. Diagram illustrating the nucleation of a bubble of true vacuum. The system is prepared 
to be in the metastable false vacuum at some early time t = —to. A bubble of size ro nucleates at rest 
at time t = 0, and subsequently expands. The semiclassical picture of an expanding bubble is valid 
for t > Tn. 



It should be noted that the process of bubble nucleation cannot itself be Lorentz in- 
variant. To illustrate this point, we may consider the limit when the bubble walls are thin 
compared to the bubble size. In this case, the trajectory of a vacuum bubble of radius r as 
a function of time t is given by [2, 3] (see Fig. 1) 



t'^ri 



(1.1) 



Eq. (1.1) is invariant under Lorentz boosts, and describes a bubble which contracts from 
infinite size to the minimum size rg, and then expands again to infinity. However, only the 
expanding part of this trajectory is relevant for describing vacum decay. Indeed, we are 
assuming that the system is prepared in the false vacuum in the far past, remaining in that 
state for a very long period of time. A tunneling event occurs at the time which we here 
denote as t = 0, forming a critical bubble, which then expands according to (1.1). Nucleation 
is a quantum process (indicated in Fig. 1 by a wavy line), and therefore we cannot say that it 
takes place exactly at the turning point hypersurface t = 0. Rather, the semiclassical picture 
of an expanding bubble is valid only after a certain time Tq has ellapsed. 



t > Tg > 0. 



(1.2) 



As we shall see, in the regime where the action of the instanton describing vacuum decay is 
large, we have Tq <^ tq. In this sense, it is still quite accurate to say that nucleation takes 
place on a t ^ const, hypersurface. 

Since only the expanding branch of (1.1) is relevant, the actual process of vacuum 
decay is not at all Lorentz invariant. In a reference frame S^ which moves at high velocity 
V = tanh (l)y with respect to the rest frame of nucleation, things look rather different (see Fig. 
2). Observers in the new frame will see a piece of the bubble appear at time t' ^ — ro sinh^-^, 
moving very fast in the boost direction. The bubble then comes to a halt at t' = 0. At 
that time the bubble wall presents a somewhat awkward shape: it consists of one hemisphere 
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Figure 2. Bubble nucleation in a reference frame S^ which moves at high velocity v = tanh^^ with 
respect to the frame of nucleation. Observers in the boosted frame will see a piece of the bubble appear 
at time t' ^ — ro sinh^^, moving very fast opposite to the boost direction. At t' = the bubble wall 
has a somewhat awkward shape: it consists of one hemisphere attached to a fuzzy interface between 
false and true vacuum, which cannot be described semiclassically. The process of formation lasts 
untile the time f ~ +ro sinh^^, when the semiclassical bubble wall closes into a full sphere. 



attached to a rather fuzzy interface between false and true vacuum, which cannot be described 
semiclassically. The nucleation process does not conclude until the time t' ^ +rosinh(/)^, 
when the semiclassical bubble wall wraps the full sphere. 

Given its impact on the kinematics of bubble formation, it is somewhat surprising that 
investigation of the frame of nucleation has been neglected for over three decades. Recently, 
however, this issue was addressed in Ref. [1]. The strategy which was adopted there was 
to consider a model detector which interacts with the bubbles, so as to reveal their state of 
motion at the time of interaction. Several plausibe scenarios were anticipated to possibly 
emerge from this study: (A) It could be that the frame of nucleation simply coincides with 
the rest frame of the detectors. In other words, each detector will see bubbles forming at rest 
in her own rest frame. (B) It is conceivable that the contracting part of the bubble history is 
not completely cut off and can be at least partially observed. (C) A third possibility is that 
the frame of nucleation is influenced by how the decaying false vacuum was set up. 

The investigation in [1] was done by using pair production of a charged scalar field (/) in 
a constant electric field as a model for bubble nucleation in (1+1) dimensions. To take care 
of option (C), the initial hypersurface where the quantum state for the field (/) is prepared 
was taken to the infinite past t ^ —oo. In this idealized situation, it was shown that the "in" 
vacuum is Lorentz invariant^, and therefore cannot determine a preferred frame of nucleation. 

To investigate (A) and (B), the model detector was chosen to be a particle of a second 
charged field i(j (see Fig. 3), interacting with (j) through the vertex 



^(0^*X + h.c.), 



(1.3) 



^This does not follow trivially from the boost invariance of the electric field, since the "in" vacuum must 
be defined in a given frame and in a specific gauge. 
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Figure 3. The "detector" particle ip can anihilate the (j) antiparticle through the vertex (1.3), pro- 
ducing a neutral particle x- Kinematically, there are at most two opportunities for the interaction to 
take place along the hyperbolic trajectory of (j). If the pair nucleates in the rest frame of the detector 
particle ?/^, as in (a), then the collision will take place in the expanding branch of the hyperbola, and 
the momentum of the decay product x will be negative p < 0. On the other hand, if the pair nucleates 
in a frame which is highly boosted with respect to the detector, then there is a good chance that the 
interaction takes place in the contracting branch of the hyperbola. This will lead to a x particle with 
positive momentum p > 0, as in (b). Therefore, a strong asymmetry in the momentum distribution of 
the products towards negative momenta can be interpreted as evidence that the detector encounters 
the expanding branch much more often than the contracting branch, consistent with option (A) in 
the introduction. 



where ^ is a coupling constant. Through this interaction, the '0 particle can anihilate the 
(j) antiparticle in the pair, producing a neutral particle x- The kinematics of this process 
is such that the (j) antiparticle has at most two chances, along its hyperbolic trajectory, of 
interacting with the detector. The reason is that the center of mass energy of the collision 
between 0* and the detector particle '0 has to be equal to the rest mass m^ of the product. 
This selects the magnitude of the momentum of (j) relative to '0, but we have two options for 
its sign. 

If the pair nucleates in the rest frame of the detector, as in the left panel of Fig. 3, then 
the collision will take place in the expanding branch of the hyperbola, and the momentum 
of the decay product x will be negative p < 0. On the other hand, if the pair nucleates in a 
frame which is highly boosted with respect to the detector, then there is a good chance that 
the interaction may take place in the contracting branch of the hyperbola. This will lead 
to a X particle with positive momentum p > 0, as in the right panel in Fig. 3. The results 
of Ref.[l] showed a strong asymmetry in the momentum distribution of the decay products, 
towards negative momenta. This was interpreted as evidence that the detector encounters 
the expanding branch much more often than the contracting branch, consistent with option 
(A). It was therefore concluded that nucleation takes place preferentially in the rest frame of 
the detector. 

The analysis of Ref. [1] was restricted to a range of parameter space where the momen- 
tum of the (j) antiparticle at the time of collision is highly relativistic. In this sense, the frame 
of nucleation was probed rather imprecisely, with a tolerance much larger than the size of 



the critical bubble. This still left some room for option (B) to be realized. 

The purpose of this paper is to sharpen the discussion given in Ref. [1], by tightening 
the precision to which the rest frame of bubble nucleation can be determined. 

We start in Section 2 by investigating the timescales which are relevant to bubble 
nucleation. These include the time Tq of quantum fuzziness after nucleation, the size of the 
critical bubble rg, the timescale Tnuc that it takes for a small quantum fluctuation in the 
false vacuum to tunnel into a critical bubble, the time to ellapsed since the false vacuum was 
prepared, and the average lifetime Tyac of the false vacuum. The hierarchy between these 
scales as a function of the instanton action Se will be clarified. 

In Section 3 we briefly review pair production by an electric field in (l+l)-dimensions. 
In Subsections 4.1 and 4.2 we consider the response of the model detector (1.3), extending the 
analysis of Ref. [1] to the range of parameter space where the momentum of the (/) antiparticle 
at the moment of collision is only mildly relativistic or non-relativistic. This will allow us to 
probe the frame of nucleation on scales smaller than tq and down to the minimum resolution 
scale At ^ Tq. Beyond that, the semiclassical picture breaks down. In Subsection 4.3 we 
use a model based on a four-point interaction, where the detector particle is not charged. 
In this case, we are able to show that even if the interaction takes place within the time 
interval At ^ Tq in the vicinity of the turning point, there is still a very strong asymmetry 
in the momentum transferred from the nucleated particles to the detector, in the direction 
of expansion after the turning point. Our conclusions are summarized in Section 5. 

2 Timescales 

In the thin wall limit, the action for a vacuum bubble is given by 



S=- f M{r)Vl-r^ dt + e j V{r) dt. (2.1) 



Here, M{r) is the mass of the domain wall of radius r, e is the difference in energy density 
between the false and the true vacuum, and V{r) is the volume inside the bubble. In D 
spacetime dimensions, we have 

M{r) ^aQ r^'^, V{r) = —^ r^-\ (2.2) 

where a is the wall tension, and ft is the surface area of the unit {D — 2)-sphere. A vacuum 
bubble has zero energy (relative to the false vacuum configuration without the bubble), and 
hence it satisfies 

Vp^ + M2 - e V{r) = 0. (2.3) 

Here the radial momentum is given by 

p^ = jMr, (2.4) 

with 7 = (1 — r^)"^/^. Up to temporal shifts, the solution of (2.3) is given by Eq. (1.1), with 

ro = ^ ^. (2.5) 

Eq. (2.5) gives the size of the critical bubble, which is also the acceleration timescale, which 
is needed for the bubble wall to become relativistic 

race ^ro. (2.6) 



To relate ro to other relevant scales, we introduce the dimensionless combination 

7 



A = 4^1 ~ 5b > 1. (2.7) 



Up to a numerical coefficient of order one, this coincides with the instanton action Se, which 
can be calculated by substituting the Euclidean version of (1.1) into the Euclidean version of 
(2.1). The semiclassical approach to tunneling is only valid when the last strong inequality 
is satisfied. This will make A a useful expansion parameter. 

Eq. (2.3) can be thought of as the classical limit of a Schrodinger equation for the wave 
function of the bubble. This was used in [2] to estimate the probability that a bubble of 
vanishing size may tunnel to the critical size. The estimate in [2] is in qualitative agreement 
with the decay rate per unit volume F ^ e~^^ which is obtained by the instanton methods 
[3]. The lifetime Tyac foi" the false vacuum is therefore exponential in A, 

Inr^acOcA. (2.8) 

As mentioned in the introduction, the semiclassical description is not adequate at the 
classical turning point r = ro, where the radial momentum pr vanishes. However, it does 
become very accurate for t ^ Tq^ where Tq can be determined from the condition [r{Tq) — 
^o] Pri^q) ^ 1- This leads to 

Tq - A-i/Vo < ro. (2.9) 

Hence, the intrinsic uncertainty in the time of nucleation is much smaller than the bubble 
size (as illustrated in Figs. 1 and 2). 

The uncertainty Tq should not be confused with the timescale Tnuc that it takes for the 
critical bubble to form out of a tiny false vacuum fluctuation. In the semiclassical picture, 
it is somewhat unclear how to characterize Tnuc- Causality suggests that this time should at 
least be as large as the size of the critical bubble. 



^nUC rsj 



ro. (2.10) 



The scale Tnuc will be characterized more precisely in the next Section, where we consider 
pair production as a model for vacuum decay in (1+1) dimensions. 
In summary, we are led to the following hierarchy of scales, 

^g < ^0 < Tnuc ^ ^0 < Tyac (2.11) 

where the strong inequalities are parametrically enhanced the larger is A ^ Se- In (2.11), 
to is the time ellapsed from the hypersurface of initial conditions to the moment when the 
bubble nucleates. This should be at least marginally larger than Tnuc We impose the last 
strong inequality on to because our primary interest is in bubbles which nucleate in isolation, 
without interference from collisions with other bubbles. Bubbles that form at to <C Tyac are 
very likely to remain isolated for a period of time which is much larger than all other scales 
involved in the problem. Finally, it should be noted that Tyac grows exponentially with A, 
whereas for all other scales the dependence is power law. In this sense, we can practically 
think of Tyac as infinite. 

In the (l+l)-dimensional example which we consider in the following Section, the electric 
field producing the pairs is external and the pairs do not interact with each other. In this 



idealized situation, the scale Tyac does not play any role, and the relevant parameters in the 
above discussion are given by 



777/ 

a = m, e^eE, A = — ;, (2.12) 

eE 

where 777. is the mass of the charged particles that nucleate, e is their charge and E is the 
electric field. 

3 Pair production 

Here we briefly review Schwinger pair production, which serves as a model for bubble nu- 
cleation in (1+1) dimensions. The main advantage is that the nucleated pairs are treated 
fully quantum mechanically. In the following Sections, we shall consider the detection of the 
nucleated pairs. Our conventions will follow those of Ref. [1]. 

Consider a charged scalar field 0, coupled to an external electric field with gauge po- 
tential 

A, = S^^A{t). (3.1) 

By spatial homogeneity, we can separate (j)(x) into Fourier modes, which satisfy the equation. 

4>k + wlcl)k = 0, (3.2) 

where 

wl ^m^ ^{k-eAf. (3.3) 

Here, 777. is the mass of the field (j) and e is its charge. 

For A{t) — —Et, we have a constant electric field £", and (3.2) becomes 

</.;' + (A + 22)0fc = O, (3.4) 

where 

777,^ 



z = VeE{t + k/eE), ^ = -^' (3-5) 

and primes indicate derivatives with respect to z. Eq. (3.4) has normalized positive frequency 
solutions 

'^' = (2jjV^''^''*^'^*["^ '""^'^]' ^^-^^ 

where 

1 + ^A 
^ = ^. (3.7) 

From the asymptotic expansion of the parabolic cylinder function at large negative z, we 
have 

-V2zy e^'^\ for - z > |zy| . (3.8) 



"^^ (2e£;)V4 
Therefore, leaving aside an irrelevant phase, we have 

<l>k^{2wk)-'/^e-'S'-"^dt' (i^_oo). (3.9) 

The "in" positive frequency mode 0/^ can be expressed as 

<t>k^ak<t>r'+pk<PT*^ (3.10) 



-i J^ Wk{t')dt' 




Figure 4. The effective potential in the Schrodinger Eq. (3.4). The scattering which determines the 
Bogolubov coefficients is overbarrier, so the reflection coefficient ^k is exponentially suppressed. 



where (j)^^{z) = {2eE)~^/^ e~^^^ Dy[{l + i)z\ is the positive frequency "out" mode, which 
behaves as the right hand side of (3.9) at late times, t -^ +oc. 

The number density of "out" particles in the "in" vacuum per unit momentum space 
interval is given by 

dn 1,^2 



dk 2ti 



m 



(3.11) 



Eq. (3.4) can be thought of as a Schrodinger equation with potential Veff = — (A + z^), and 
the determination of the Bogolubov coefficients ak and ^k amounts to solving a scattering 
problem in this potential (see Fig. 4). Note that k enters Eq. (3.4) only through the combi- 
nation k + eEt. This means that the scattering near z = occurs approximately at the time 
t ^ —k/{eE), but it is otherwise independent of k. Eq. (3.11) can then be rewritten as 



r = dn/dt = eE\(3k\'^/{27T) 



(3.12) 



where the coefficient (3^ = e^^^ can be easily calculated from a standard linear relation 
between the parabolic cylinder functions appearing in 0/^ and 0^^* (see e.g. [1] and references 
therein). This leads to Schwinger rate 



eE 

2^' 



-nX 



(3.13) 



The exponent in (3.13) coincides with the Euclidean action Se = ttA, in agreement with the 
instanton approach, which is valid for large Se- 

The asymptotic expansion (3.8) for (1)^, and the analogous one for 0^^^, are valid for 
1^1 » \iy\. Hence the timescale of scattering off the potential in Fig. 4, where the mixing of 
positive and negative frequency modes occurs, is at most of order 



^mixing n^ I'yCrj) 



-1/2 ^ ;^l/2 



^'O- 



(3.14) 



On the other hand, this seems to be a rather crude upper bound. Mode mixing wih be at its 
peak when the non-adiabaticity parameter 



fk{t) 



Wk 



wt 



(3.15) 



is at its maximum. The function fk is symmetric around the time t = tk = k/eE (where it 
vanishes) and has two peaks which are at a distance At ^ rg away from t^- Hence, we can 
estimate 

rmixing ^ r^. (3.16) 

Identifying this with the semiclassical timescale Tnuc foi" nucleation of a pair out of a vacuum 
fluctuation suggests that 

Tnuc ^ rQ. (3.17) 

This is consistent with the estimate Tnuc ^ '^o which we mentioned in the previous subsection. 
Stronger evidence for (3.17) can be found by considering the case of an electric field 
E(t) which is turned on and off on a timescale to: 



m 



E 



cos\i\t/to) 



(3.18) 



The corresponding gauge potential is given by A{t) — — E'tg tanh(t/to), and the mode equa- 
tion (3.2-3.3), can in this case be solved in terms of hypergeometric functions . The Bogoli- 
ubov coefficient is given by (see e.g. [7] and references therein) 



Wk? 



cosh 



T,J{eEtlY-\ +smh2[^(«;+-u;_)] 



sinh(7rtoit;+) smh{TitoW- 



(3.19) 



where w± — \Jw? + (A: ^ eEt^^ . For the mode with conserved momentum A: = 0, the 
physical momentum k^hys — k — eA{i) vanishes precisely at the time t = 0, which therefore 
corresponds to the turning point. For this mode, Wj^ — W- and (3.19) simplifies to 



\/3k- 



cosh 


2 


TT^ieEtir - I 




sinh^ 


TT 


to^m'^ + (eEto) 


2 



(3.20) 



This expression determines, through Eq. (3.12), the rate of pair production at the time when 
the electric field is at its maximum E{t — {)) — E. Further simplification occurs when the 
timescale to during which the electric field is switched on is bigger than the scale of quantum 
fuzziness t^^ Tq ^ \~^'^r(). In this case, we have 

eEtl > A^/^>1. {Tq<t^). (3.21) 

Now, we can distinguish two cases. For tg <C tq, we have eEt^ <C m, which leads to 

|/3fc=o|' ^ e-2'^'"*o, (Tg < to « ro), (3.22) 

while for to ^ ro, we have eEt^ » tti, and 

|/3fc=oP«e-^\ (to>ro). (3.23) 



In the first case, Eq. (3.22), the pair production rate near t = is independent of the electric 
field. Clearly, this is unrelated to the Schwinger process, and corresponds to particle creation 
by the non-adiabaticity of the switching process. On the other hand, we see from Eq. (3.23) 
that if the electric field changes on timescales which are larger than the radius of the critical 
bubble, then we obtain the adiabatic Schwinger pair production rate. This strongly suggests 
that the timescale of nucleation of a pair from a small vacuum fluctuation satisfies Tnuc ^ ^o- 
In combination with the bound Tnuc ^ ^o which is suggested by causality, we are again lead to 
the estimate (3.17). In conclusion, out of the scales in (2.11), two of them are of comparable 

size ro ^ Tnuc- 

4 Charged detector (3-point interaction) 

Let us start by reviewing the model detector which was introduced in Ref. [1]. This involves 
a charged '0 particle interacting with the antiparticle of a nucleated pair of the charged field 
0, through the vertex 

^((/)^*X + h.c.), (4.1) 

In the gauge (3.1), the canonical momentum of charged particles is conserved. Denoting by 
—k the momentum of the (/) antiparticle, by q the momentum of the '0 particle, and by p the 
momentum of the x particle, we have 

p = q-k. (4.2) 

The physical momentum of charged particles changes with time, due to the acceleration 
caused by the electric field, 

Qphys = ^ + eEt, kphys = A: + eEt. (4.3) 

From (4.2) and (4.3), we also have p = Qphys ~ kphys- 

Note that the tp particle detector is accelerated by the electric field. At early times it is 
highly relativistic with negative velocity, and at late times it is highly relativistic with positive 
velocity, following a hyperbolic trajectory. As a consequence, the momentum distribution of 
the resulting x particles is actually Lorentz invariant [1]: 

dTV = C^, (4.4) 

where AL is the number of x particles, p is their spatial momentum. 



w^ = Jp'^ + m^ (4.5) 



and C is a constant. The form of (4.4) can be understood as follows. The interaction of tp 
with nucleated pairs has a constant probability of occurring per unit proper time interval dr 
along the trajectory of the tp particle. This trajectory can also be thought of as an orbit of 
the Lorentz boosts, with parameter (py = ctipr^ where a^ = eE/m^ is the proper acceleration 
of the detector particle. If we consider two collisions which look identical in the detector's 
rest frame, separated by a time delay dr, the corresponding momenta p will be related by 
a boost with parameter dcpy = (dp/w^) = a^dr. Hence, Eq. (4.4) simply expresses the fact 
that the probability of producing a x particle per unit proper time is constant. An immediate 
consequence is that the distribution (4.4) cannot inform us about the frame of nucleation. 
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This difficulty is avoided by switching on the interaction (1.3) for a short period of time 
T. This can be implemented by substituting the coupling g in (4.1) by a time dependent 
coupling with a Gaussian proffie [1]: 



g ^ g{t) = g e"*'/'^' 



If we choose 



then, during the time interval |t| < T", we have 



kphysl ^ eET <^m^ <m^- m^. 



(4.6) 
(4.7) 

(4.8) 



For any time interval T, we may choose the mass m^ of the detector to be heavy enough, 
so that the above strong inequality holds, and in this case, the detector is essentially at rest 
within the time interval T before the interaction. The last inequality in (4.8) indicates that 
we are choosing the mass m^ of the decay product to be larger than m^ + tti^^. Otherwise 
the detection of (j) antiparticles would be forbidden by energy conservation (at least in the 
absence of an electric field and of any time dependence in the coupling g). Within this setup, 
a significant asymmetry in the momentum distribution for the decay products x towards 
negative values of p will indicate that the detector is more likely to be hit from the right, 
indicating that the interaction is more likely to happen in the "expanding" branch of the 
hyperbola, when the (j) antiparticle is moving to the left. On the other hand, a more symmetric 
distribution would indicate that there is a significant chance of detecting the anti-particle 
also in the "contracting" branch (see Fig. 3) . 

The momentum distribution of x particles for an interaction of finite duration T was 
calculated in [1], and it is given by 



dN. 



where 



^ = ^l^(p.. = o,r)|^ 



/oo 
dt g{t) Vg=o(t) </'-p(t) x;(t)- 
-OO 



(4.9) 



(4.10) 



Here, the mode function 0/^ is given by (3.6), and similarly for '^g (with the mass tti^^ re- 
placed with m^). The mode function for the neutral particle is simply given by Xp — 
(2i(;^)~-^/^e~^^^^. It follows from Eq. (4.8) that the ijj particle stays non-relativistic through- 
out the interaction, so we can approximate '0o ^ {2m^)~^''^e~^'^e~'^^'^'^ . Here, particle 
production of the '0 field is also neglected. Within this approximation, the amplitude (4.10) 
can be calculated explicitly [1], and we have 



dN^ Ct 
^ ^ — exp 
dp w^ 



\4 + a^)eE 



\DAZ)Y 



where 
and 



z/ = -(l + zA)/2 
W = w^ — m^. 



(4.11) 

(4.12) 
(4.13) 
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f<^T, 




p<0 

Figure 5. The left panel shows the momentum distribution of x particles, Eq. (4.11), for m^/y = 
5 (e£;)i/2 (corresponding to A = 25), with m^ = 110 {eEy^'^, m^ = 100 (eE)^/'^ and a = eET'^ = 
A^/^. The width Ap of the distribution can be estimated from (4.8). During the time interval of 
duration T when the interaction is switched on, indicated by shading in the right panel, the detector 
particle is accelerating and its momentum changes by the diinount Aqp^ys ~ eET. This leads to a 
similar spread Ap ^ eET in the momentum distribution of the products. The distribution is clearly 
biased towards negative momenta p < 0. This means that the detector ip never encounters the (p 
antiparticle in the contractig branch of the hyperbola (represented by the dotted line). 



Note that if energy conservation were exact, then W would be the energy of the charged 
antiparticle (f). In (4.11) we have also introduced the parameter 



a = eET'^, 



(4.14) 



and the variable 



Z = 



2 aW + 2ip 
JeE 



V4T 



a^ 



Finally, the momentum independent prefactor in (4.11) is given by 



Ct 



9^ e 



-37rA/4 



a 



2m^(2e£;)3/2^4ip^ 



3Atan-i(2/Q;) 



(4.15) 



(4.16) 



A plot of the distribution (4.9) is given in the left panel of Fig. 5, for m^^ = 5 {eE^^^ 

(corresponding to A = 25), with m^ = 110 [eE^I^ , m^ = 100 (eE)^/'^ and a = A^/^. It is 

clear that the momentum distribution is highly biased towards negative momenta, p < 0. 

Denoting by tc the time of collision minus the time at which nucleation occurs, the 

physical momentum of the (f) antiparticle at the time of collision is given by p = —eEtc, and 

we have 

P ^ _^ 

m^ ro' 

Consider, for illustration, the distribution in Fig 5. This is peaked at p^ ^ m?,^ which 
corresponds to a time of collision \tc\ ^ vq (see the right panel of Fig. 5). For the chosen 
values of the masses tti^^, m^ and tti^, both signs oi p and tc would be kinematically allowed, 
but we find that detection only occurs in the expanding branch of the hyperbola. From this 
we can conclude that the boost parameter (^^ which relates the frame of nucleation and the 
frame of the detector (see Fig. 2) is bounded by 



(4.17) 



sinh|(?^„| < 






(4.18) 
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This sets the accuracy to which the frame of nucleation is determined by this particular 
interaction. In the example represented in Fig. 5, we have tc ^ tq, which leads to \(f)y\ < 1. 
This is not too accurate, but we shall see that the precision can be increased parametrically 
in A, by choosing m^ — m^ to be sufficiently close to the rest mass of the antiparticle (/). 

Our goal here is to examine the regime where both the time of collision tc and the time 
interval T during which the interaction is switched on are marginally bigger than, but very 



close to the scale r^. 



t^r^T r^ (few times) Tq < tq. (4.19) 



In this way, we can probe the time of nucleation to the finest possible resolution consistent 
with the uncertainty principle. To realize this situation, we need to explore the regime of 
very large A. The reason is that the hierarchy between the quantum scale and the size of 
the critical bubble, Tq ^ A~^/^ro, is only appreciable for, say, A > 10^. In this parametric 
range, the order of the parabolic cylinder function in (4.11) is very large \u\ > 10^. As a 
result, a numerical study becomes quite cumbersome, and an analytic approach will prove 
more illuminating. 

First, we note that from (4.8), and (4.13) we have 

W > m^, (4.20) 

and using (4.15), we have 

|Z|>^>aV2»i. (4.21) 



The parabolic cylinder functions of large order |z/| ^ A » 1 and argument Z > A^/^ can be 
estimated by using an integral representation of D^{Z)^ and then evaluating this integral by 
using the method of steepest descents. This is done in Appendix A. In this approximation, 
the distribution (4.11) is given by 



dNy ^ — -^ exp — 2- 

^ 8eE m^ P ^ 



{eETf 



^. (4.22) 



Here, we have introduced 



{W^-ml)^/\ (4.23) 



Note from Eq. (4.13) that W = w^ — m^. Since, from (4.8) '0 is assumed to be non- 
relativistic, W is the difference in energy between the x ^^^ the tp particle. Hence, if energy 
is approximately conserved, we have W ^ w^f) = {p'^ + mV)^/'^^ and so we expect |p| ^ 
P, in agreement with (4.22). Even though P is strictly speaking p dependent, it will be 
approximately constant for p <C m^. The mass of the x particle can be taken as large as 
desired, so we can always think of 



as a constant parameter in the momentum range of interest. In deriving (4.22) we have also 
assumed, according to Eqs. (A. 18) and (A. 19), that 

P > A-^/'m^, (4.25) 
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and 

P > a-^m^, (4.26) 

with 

a = eET'^ > 1. (4.27) 

Eq. (4.22) is the main result of this Section. 

The momentum distribution (4.22) is peaked at negative values 

p^-P <0, (4.28) 

with a width 

Ap - eET. (4.29) 

The origin of the spread in momentum in Eq. (4.29) is simple to understand: it corresponds 
to the variation in the momentum of a charged particle (the detector) due to the force exerted 
by the electric field E during the time interval T. 

We are particularly interested in small values of the momentum, for which (4.25) is 
marginally satisfied, 

\P\ - A-^/^m^. (4.30) 

In this case, the condition (4.26) requires a » A^/^, or 

T > A-i/Vo - Tq. (4.31) 

On the other hand, the asymmetry in the momentum distribution will only be appreciable 
provided that \P\ ^ A"-'^/^??!^^ > Ap ^ eET. This requires 

T < A-V3^o - Tg. (4.32) 

Hence, we can have a well defined peak at p ^ — A"^/^??!^^ by choosing the mass difference 
m^ — m^ in (4.24) so that P satisfies (4.30), and by using T ^ Tq. Combining with (4.17) 
and (4.18), we find that the boost parameter (/)y which relates the rest frame of the detector 
with the rest frame of bubble nucleation can be parametrically bounded with precision 

corresponding to a time of collision tc ^ A~^/^ro ^ Tq. Hence, the larger the Euclidean action 
Se^ the sharper is the precision to which both rest frames coincide. 

To understand the origin of (4.33), we note that energy conservation at the time of 
collision becomes imprecise in inverse proportion to the time of interaction At = T. For 
non-relativistic particles, we have 

^ > {At)-\ (4.34) 

2m^ 

Combining (4.29) and (4.34), the minimum uncertainty in the velocity of the (j) particle is 
given by 

Av > A-i/3, (4.35) 

which is saturated by At ^ Tq. This uncertainty in the velocity of a charged particle embedded 
in an electric field cannot be reduced. 
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Note also that if we use a 3-point interaction such as (4.1), and we take the neutral 
field X ^s the detector, then the only product of the interaction is a charged particle '0 
which accelerates uniformly. Because of that, if the interaction is not switched on and off, 
the resulting distribution of physical momenta for '0 is Lorentz invariant [1], and shows no 
asymmetry in momentum. Hence, we cannot tell whether the detector has been hit from the 
right or from the left. The only way to check for this asymmetry would be by considering 
an interaction of finite duration T. But in that case we encounter the same situation which 
we discussed above: the momentum at the time of collision is uncertain by the amount given 
in Eq. (4.29), and the precision in the determination of the frame of nucleation is limited by 
(4.35). 

5 Neutral detector (4-point interaction) 

To avoid some of the difficulties associated with the time variation in the momentum of 
charged particles, in this Section we consider a 4-point interaction 

^(XiX20^* + XiX20V), (5.1) 

where both the detector particle xi ^^d one of the products of interaction X2 are neutral 
with respect to the electric field. The main advantage of this setup is that we will not need to 
switch on and off the interaction in order to find an asymmetry in the momentum distribution 
of the products. 

Even though xi ^tnd X2 ctre neutral, we take them to be complex, so that they carry a 
conserved global charge. The reason is that we want to detect (/) particles, but not (j) anti- 
particles. With the interaction (5.1) this is achieved by using an incoming xi particle as our 
detector (as opposed to a xi antiparticle) . For simplicity, we choose 

rui = 1712 = '^x- (5-2) 

This automatically suppresses the direct decay of the xi particle into the other three (see 
Fig. 6). 

The expectation value in the number of X2 particles can be calculated by using the same 
techniques which were used in Ref. [1] for the case of the 3-point interaction. This leads to 
the expression 

where the amplitude is given by 

/oo 
dt' <Pl(t') r,{t') Xi,o(i') X*2,p{t'\ (5-4) 

-OO 

with 

k — q^ p. (5.5) 

Without loss of generality, we have taken the initial xi particle to be at rest, p' — 0, while 
X2 has momentum p — p. For the neutral particles Xi^ ^he mode functions in the amplitude 
(5.4) are given by 

Xi,o = {2m^)-^/^e-'^-\ X2,p = {2w2)-^'''e-'''^K (5.6) 
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(a) 



(b) 



Figure 6. The four-point interaction. Panel (a) represents detection of a ^ particle by interaction 
with XI7 while panel {b) represents direct decay of xi- The latter process is suppressed for mi < 
^2 + ^0 + ^^, which is guaranteed for mi = m2 = m^. 



where 



W2 



{p^ + rr? 



a/2 



(5.7) 



For the charged particles, the mode functions '^g and (j)k are given by Eq. (3.6), where the 
corresponding values of the masses and momenta should be used in the definitions of z and 
AinEq. (3.5). 

5.1 Kinematics 

Since we have chosen mi — 7712, the neutral sector can be thought of as a detector of mass 
m^ which is being hit by a charged particle (j). During the interaction, part of the momentum 
of this charged particle is transferred to the detector, and part of it to a charged particle 
'^, whose mass is m^ > rricf). It is convenient to introduce the invariant momentum transfer 
(which is space-like). 



The physical momenta of the charged ^ and (j) particles are given, respectively, by 

(Iphys — ^ + eht^ f^phys = A: + eht. 



(5. 



(5.9) 



Because of the time dependence in (5.9) we cannot properly speak of energy conservation at 
the collision. Rather, as mentioned around Eq. (4.35), there will be some uncertainty in the 
conservation of energy which depends on the time duration At of the interaction AEAt ^ 1. 
In the non-relativistic regime, the uncertainty in the momenta is minimized by choosing 



At r^ Tq, in which case 



^kphys = Mphys ^ A ^/^m^. 
Within such uncertainty, the physical momentum 

Kp = i^physy^c) 



(5.10) 
(5.11) 
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of the (j) particle at the time of colhsion can be uniquely reconstructed from the value of the 
spatial momentum p — p oi the X2 particle. Indeed, the approximate energy conservation 
equation 



I ml + kl + m^^ ^m\ + {k^ - pf + ^m\ + p^, (5.12) 

can be written as 



ml + fc2 ^ ^ml + (fc, -p)2 + ^— (5.13) 



U_ 



It is then straightforward to show that 



fcp « ^ + sign(p)-V — "-, (for kp > A-V3^^), (5.14) 

where we have introduced 

M^ = ml-ml-hU. (5.15) 

It can be checked that the double sign in the square root should be taken as indicated in 
(5.14), so that Eq. (5.13) holds before we square it. Thus, we find that there is a single 
value of kp for each value of p, and a measurement of p amounts to a measurement of the 
momentum of the charged particle at the time of collision. It is also clear from (5.14) that 

sign(A:p) = sign(p). (5.16) 

This is convenient because the sign of p tells us whether our detector has been hit from the 
left or from the right. 

5.2 Asymmetry 

Let us now look at the distribution (5.3), to check for an asymmetry between positive and 
negative p. First of all, we note that '0 particles are charged and they can also be pair- 
produced. However, we are only interested in detecting (j) particles, without contamination 
from ^ particles. Hence, we shall take 

A^ - A^ > 1, (5.17) 

where we have introduced 

2 

7T7 . 

A. = ^. (5.18) 

Eq. (5.17) ensures that the rate of production of '0 particles is much smaller than the rate of 
production of (j) particles. The amplitude (5.4) is calculated in Appendix B.l. Introducing 
(B.15) in (5.3), and using Eqs. (B.16-B.19), we have 



dtdp ~ 16^ V ^"""^^ / ^2m^U 



W-ia^^ia^ I i 



—) 

2eEJ 



\ for{^;»: ,5.,) 



Here, 

a± = ^^^. (5.20) 

In deriving (5.19), we have taken into account that \A2\ is independent of q. From (5.9), 
this can be understood in terms of the invariance of the interaction probability under time 
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translations. For given Qphys the value of q is different depending on the time at which the 
interaction takes place, so we have dN2/dtdp = (eE)dN2/dqdp. 

The most remarkable feature of Eq. (5.19) is that, together with (5.17), it shows a strong 
asymmetry in the distribution of p between positive and negative momenta: 



dN2{-\p\) 



If we combine (5.21) with (5.16), this implies a strong asymmetry in the number of times 
the collision happens with positive or with negative kp. The smaller is kp, the closer the 
collision is to the turning point of the particle's hyperbolic trajectory, and the better we 
test for the frame of nucleation. In practice, the fundamental uncertainty Akp ^ A"^/^??!^^ 
in the momentum of a non-relativistic charged particle will limit the precision to which we 
can determine the relative velocity between the frame of nucleation and the frame of the 
detector, so that Av ^ A~^/^. Nonetheless, (5.21) is valid beyond the semiclassical regime, 
and the distribution of the momentum transfer can be very strongly biased in the direction 
of the expanding branch, even if the interaction takes place very close to the turning point. 
Let us now show that there is a parametric range where the asymmetry (5.21) is quite 
significant, while the collision cannot be described semiclassically. First, we note that the 
condition (5.17) for a sizable asymmetry also implies that there is a minimum value of kp, 
which we may call /cq, for the collision to be kinematically allowed. To find out the value of 
fco, we can calculate the square of the total 2-momentum before and after the collision. The 
minimum value of kp is such that the decay products have no relative speed, so that the total 
energy is the least possible. Then, in the center of mass frame, the total 2-momentum P^ of 
the products of collision only has temporal component, and we have 

P^P^^(m^ + m^)2, (5.22) 

where the equality is only approximate because of the fundamental uncertainty in the energy 
of charged particles. Before collision, the total 2-momentum in the frame of the xi particle 
is given by P^ ^ (tti^ + a/tti'^ + /cq, /cq ) . Squaring and equating to (5.22) we have 



7 2^/2 2> 



^^m^^m^-m^ 



m^ 4771^ 



(5.23) 



It is clear that ko will be minimized by taking the detector particle to be much heavier than 
any of the other particles 

m^ » m^, 171(1). (5.24) 

In this case, we have 

fc2 ^ml-ml> A- V2, (5.25) 

where the approximate inequality is due to (5.17). From (5.14), this corresponds to momen- 
tum transfer of the order 

pl^kl^ml- m|. (5.26) 

— 1/3 

On the other hand, the uncertainty in the value of ko is of order A A; ^ A , tti^^. 
We conclude that by choosing m^ to be in the range 

A^'ml < {ml - ml) < \-^^'^m% (5.27) 
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and for momentum transfer of the order p^ ^ Pq, the semiclassical approximation is not vahd, 
while the asymmetry (5.21) in the distribution oip is still quite sizable. Hence, the asymmetry 
persists for interactions which occur very close to the turning point of the semiclassical 
trajectory. 

6 Summary and conclusions 

We have studied the frame of nucleation of charged particles produced by the Schwinger 
process, by using two different model detectors. 

The first model is based on the 3-point interaction (1.3) which was introduced in Ref. 
[1]. In this case, the detector is a charged particle ^/^ which anihilates with the (j) antiparticle of 
a nucleated pair, producing a neutral particle x- Semiclassically, the trajectory of a nucleated 
pair is hyperbolic, with incoming and outgoing branches separated by a turning point. In Ref. 
[1] it was shown that in the detector's frame, nucleated particles are preferentially encountered 
after the turning point, with a strong asymmetry between the number of incoming and 
outgoing particles. The analysis of [1] concentrated on the detection of particles which 
are fairly distant from the classical turning point, moving already at ultrarelativistic speed 
1 2; I ^ 1. In this sense, the relative velocity between the frame of nucleation and the frame of 
the detector was not very well constrained. Here, we have refined this analysis, extending it 
to the regime where the (f) particle is non-relativistic at the time of collision. We have shown 

— 1/3 

that the asymmetry persists down to velocities of order Av ^ S^ <C 1. Here, Av is the 
minimum uncertainty in the velocity of a non-relativistic charged particle in the presence of 
an electric field. This represents a fundamental limitation in the determination of the relative 
speed, which cannot be overcome. 

The second model we have introduced is based on the 4-point interaction (5.1). Here, 
the detector is a neutral particle xii which upon interaction with the charged particle (j) 
produces a second neutral particle X2 and the charged particle '0, 

XI + 0^X2 + ^. (6.1) 

We choose the masses of xi and X2 to be the same, so we can think of the neutral sector as a 
detector of mass m^ which is being hit by a charged (j) particle. In this case, we find a strong 
asymmetry in the momentum p transferred from the charged to the neutral sector, pointing 
in the direction of expansion after the turning point [see Eq. (5.21)]. Like in the case of the 3- 
point interaction, this means that throughout the semiclassical regime, the detector is almost 
exclusively excited by the outgoing (i.e. expanding) branch of the hyperbolic trajectory, and 
hardly ever by the incoming (i.e. contracting) branch. Moreover, the strong asymmetry 
persists even for values of p corresponding to interactions which take place near the turning 
point, where the semiclassical approximation breaks down. We conclude that the rest frame 
of nucleation is very sharply defined and strongly correlated with the frame of the detector. 
There are several directions in which our analysis could be extended. First of all, we 
have restricted attention to the case where the electric field is switched on for an infinitely 
long time, and the charged field (j) is in its Lorentz invariant "in" vacuum. More realistically, 
the electric field decays by pair production, so its time duration should be finite. In this 
case, we expect that the frame of nucleation can also be influenced by the frame where the 
electric field is turned on at some finite initial time. We expect that initial conditions will 
not play a role provided that the detector's proper time r elapsed since the hypersurface of 
initial conditions is much larger than the time Tnuc that it takes for a pair to be produced 
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from a vacuum fluctuation. In Section 3 we have estimated this time to be of the order of 
the size of the instanton tq. We therefore expect our results to be vahd for r ::^ tq. A more 
quantitative investigation of this case is left for further research. Another direction which 
seems worth pursuing is the case of pair production by an electric fleld in de Sitter space [8] . 
We expect to come back to this issue in a future publication. 
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A Asymptotic expansion of Dj^{Z) 

The parabolic cylinder functions of large order \v\ ^ \^ 1 and argument Z > A^/^ can be 
estimated by using an integral representation of Djy{Z)^ and then evaluating this integral by 
using the method of steepest descents. 
Following [9], we write 

D,{Z) = ^S^±J^e-''''Z I ef^^'^^Uv, (A.l) 

where the contour C is depicted in Fig. 7, and 

f{v; Z) = Z'^^v - v^/2) -{l + v) \n{vZ). (A.2) 

The saddle points are determined by f'{vi : Z) = 0, where a prime indicates derivative with 
respect to v. This equation is solved by 

2^. = 1 + |, (A.3) 

Here we have introduced 

S{Z) = ±[Z^-4.{l + v)Y'^. (A.4) 

In principle we should consider both possible values of the square root. However, from the 
location of the saddle points relative to the overall geometry of the contour, we shall see that 
the only relevant saddle point which is the one closest to the branch cut oi f{v; Z) (see Fig. 
7). From Eqs. (4.15) and (A.4), we have 

S' - -^(1^^ [^"P" - 4(P' + <) + ^i^PW] - 2, (A.5) 

where 
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Figure 7. The contour of integration C in Eq. (A.l). The saddle points vf^ given by Eq. (A. 15) are 
also represented. The root that is closest to C is v^, and will give the dominant contribution provided 
that Eq. (A. 10) is satisfied. In this case, the saddle point approximation is obtained by deforming C 
so that it goes through the point vj . 



By performing the Gaussian integral around the saddle point, we have 



^f^^■^^U^ 



c 



27r 



1/2 



where 



r: 



-2Z^ 



s 



z + s 

is the second derivative of / with respect to v at the saddle point, and 



(A.7) 
(A.8) 
(A.9) 



The saddle point approximation is valid when the distance \5v\ = \S/Z\ between the two 
saddle points is large enough, so that the Gaussian integrals around both saddle points do 
not have much overlap \fg5v^\ S> 1. This requires 



\S^\ > \s + z\. 



Substituting in (A.l), we have 



Di,{Z) ^ -ir(l + v) e 



(l+u)/2 



z + s 



-i^+i^) pZS/4 



V2^' 



from which we find 



\DAZ)? 



62 



|-Aarg(Z+^) 



,ZS/2 



(A.IO) 



(A.ll) 



(A.12) 



2cosh(7rA/2) 
To proceed, we shall assume that the difference m^ — m^ is large enough, so that 

aP > m^. (A.13) 
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We also assume that \p\ < P and a ^ 1. From (A. 5), this leads to 



^p\ 2^pW_^eE 
ieE ] aP^ 2P^ 



To lowest order, from (A. 3) we have the two solutions 

P 2ip ieE 






(A.14) 



2.±^l±f^+'^^ 



2ipP\ 
W ' aP 2PW aW^J 



(A.15) 



It follows that the saddle point which is the closest to the branch cut of the function f{v;Z) 
given in (A. 2) corresponds to the lower sign. It is then straightforward to check that 



- Aarg(Z + S) + ^ZS/2] « -4^ ^T'l' 
Substituting (A. 12) in (4.11), and using (A. 16), we obtain 



-ttA 



dN^ 



8eE m^ P 



exp 



, (p + P)^ 1 dp 
'{eET)^\w^' 



(A.16) 



(A.17) 



Finally, we should summarize the conditions for the validity of (A.17). Using (4.15) and 
(A.14), Eq. (A. 10) requires P^ » eE W. This is satisfied provided that 



P > A-V3 



m^. 



(A.18) 



For smaller values of P, the two saddle points are too close together and the Gaussian 
approximation to the integral representation (A.l) is not valid. Also, since we have used 
(A. 13) in order to expand S and vf^ we must require 



P » a ^rricf), 



(A.19) 



with a = eET'^ > 1. 



B Transition amplitude for the neutral detector 

B.l Probing the "in" vacuum state of the (/) field 

Here we shall calculate the amplitude (5.4) for the 4-point interaction, starting out in the 
"in" vacuum state for the (f) field. For generality, and as a check of Lorentz invariance, here 
we consider arbitrary momenta p' and p for the neutral particles xi ^tnd X2- Then we have 



M{k,q,p]p'\ 



where we have defined 
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e '^- 



{2eEy/^(2upy/^(2ujp^y/^ 



z^ = VeE t + 



eE 



l^cfy = 



1 + iAw 



A^ = 



m^ 



eE' 



(B.l) 
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., = ^4E(t+-|), ., = -i^, A, = §. (B.2) 

Using the following integral representation of the parabolic cylinder functions: 

D^[-V2e^'z] = ^ -e-t^' / dx x^* e^^^^-% (B.3) 

the amplitude becomes 



A2{k,q,p; p') 



^¥+n/^^?/-2^(^'+^') 



(2e£;)V2(2c^^~)i/2(2a;^~,)V2 r(-^^)r(-^^) 

X / dx x""^^ e'^^'^'y^ I dyy'^^e^'^^'y^'''^' 
Jo Jo 

/oo 

^^ ^-ieEt'^-i{uj~,-ujp+k+q-V2eEx-V2eEy)t m ^\ 

-oo 

After completing the square in the exponent of the t-integral and performing the Gaussian 
integration, we have 

roo /♦oo 

X dx x"'^ e-^^+^ / dy /^ e^(^-'^-)^ , (B.5) 

Jo ^0 

where we defined 

^^^ — vn — ' ''-^ — vm — • ^""-'^ 

Note that the product of /i+ and /x_ is related to the momentum transferred to the neutral 
sector at the collision 

Defining 

u — —fJi-^ /3 = i//+, (B.8) 

the y integral in Eq. (B.5) is well defined for 

< arg(ix) < TT, (B.9) 

and in that case it is given in terms of Euler's Gamma function. The x, ^-integrals in Eq. (B.5) 
then become: 

/•oo /»oo /»oo 

/ dx x^'l e-^^ / dy y^^ e^(^+^)^ = e-'^^^V{-v^) / dx x^'l (x + uf^ e-^^.(B.10) 
Jo Jo Jo 

Likewise, the existence of the integral in the right hand side of Eq. (B.IO) requires 

- ^ < arg(/3) < ^. (B.ll) 
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From Eq. 3.383.4 in Ref. [10] we have 

dx x^i {x + uY* e"^^ = iS-^-'^-u-'^+'^-Ti-u^,) e^ W-ia+,ia.{M- (B.12) 



f 

Jo 



/o 
Here, we have introduced 

a± = (B.13J 

Eq. (B.12) is valid for real positive values of the parameters /3 > and u > 0^ but it can be 
analytically continued to the domain 

3?/3>0, \d.Tg{u)\<7v. (B.14) 

If u were real and negative, the contour of integration would step over a singularity at x = —u, 
so the integral has a branch cut at i/ < 0. On the other hand, the domain (B.14) includes 
the parameter region of our interest, which is given by (B.9) and (B.ll). 
Using Eqs. (B.5), (B.IO) and (B.12), we have 



\A2{p;p')\'=9'^^^, ^^ ^^\W-^a^,,_{H\^ (B.15) 

Now, using (B.ll) and (B.9) we can disambiguate the phases of /3 and u as follows. 



and 



^=i..i(:-fi), for{;;;^2: m,) 



I"-! if). '-{::; 2: (B17) 



From (B.6) and (B.7), it is clear that /x+ and ^_ are the null components of the momentum 
transfer 2-vector, and therefore the signs sign(/i±) and the product (//+//_) are Lorentz 
invariant. Hence, Eq. (B.15) for the modulus squared of the amplitude is written in a 
manifestly Lorentz invariant form. 

In the frame where the initial particle is at rest, we have p' — {)^ p — p. Moreover, in 
the case when the masses of the neutral partilces are degenerate, mi — 1712 = m^, we have 

sign(A^+) = -sign(^_) = -sign(p), (B.18) 

while from (B.7) we have 

Im+M-I = ^>0. (B.19) 

Introducing (B.15) in (5.3), and using Eqs. (B.16-B.19), we derive Eq. (5.19) in the main 
text. 

B.2 Probing a one particle state for the (/) field 

As a consistency check, and an alternative way of deriving the momentum discontinuity in 
Eq. (5.19), let us now "calibrate" our 4-point interaction detector by introducing an actual (j) 



24 



particle of momentum k in the initial state, in addition to the neutral xi particle. Changing 
(j)k{z) to (f)l.{z) in Eq. (5.4), the new amplitude is given by 

JO Jo 



/o 



J —{ 



^^ ^-i{ujp/-u;p-k+q+V2eEx-V2eEy)t m 2]_^ 



Using q + p — k -\- p' — 0, the delta function which results from time integration can be 
written as 

/oo o 

dt e-ii^p^-^p-k+q+V2^x-V2^y)t ^ ^,^ S{ti+ +x-y). (B.22) 

For /x+ > we do the y integral, and after using Eq. (B.12), we obtain: 

while for /x+ < 0, we must do the x integral first, and we obtain 

4(e£;)^^p^^/ |r(-z/*)|^ Im+m-I 

Note, however, that Wy,^^{z) = W/^-jy{z)^ and so the factor containing the Whittaker function 
is common in Eqs. (B.23) and (B.24). Finally, setting p — p and p' — 0, we have 



1^2 1 
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f°^^p<o: ^^-'^^ 



The dominant term in (B.25) is symmetric in p, and corresponds to the fact that the initial 
(j) particle is in hyperbolic motion. For any given |p|, its trajectory has two chances for 
interacting with the detector. In one case the detector is hit from the right, and in the other 
case from the left, which gives rise to the two possible signs of p with equal probability. 
The subleading term in (B.25), on the other hand, has the same momentum asymmetry as 
the distribution (5.19). This is not surprising, since aside from the original (j) particle of 
momentum A:, we will have a bath of pairs of (j) and ^ particles which are created by the 
electric field. In particular, there will be a pair of (j) particles of momentum k which will be 
created, with probability proportional to e~^^^. Upon interaction, the created (f) particle will 
deposit positive momentum transfer p > on the detector. Likewise, a pair of '0 particles of 
momentum q can be created with probability proportional to e~^^^, and upon interaction, 
the '0 antiparticle will deposit negative momentum transfer p < on the detector. 
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